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Abstract 

The propagation of electromagnetic waves in an unmagnetized weakly inhomogeneous cold 
plasma is examined. We show that the inhomogeneity induces a gauge connection term in wave 
equation, which gives rise to Berry effects in the dynamics of polarized rays in the post geometric 
optics approximation. The polarization plane of a plane polarized ray rotates as a result of the 
geometric Berry phase, which is the Rytov rotation. Also, the Berry curvature causes the optical 
Hall effect, according to which, rays of left/right circular polarization deflect oppositely to produce 
a spin current directed across the direction of propagation. 
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I. INTRODUCTION 



It has often been said that ninety nine percent of the matter in the universe is in plasma 
state fl. Because real plasmas have density inhomogeneity, the study of electromagnetic 
wave propagation in inhomogeneous plasmas has become important [2h4|. In this paper 
we examine the propagation of electromagnetic waves in an unmagnetized plasma via the 
post geometric optics approximation {5, 6|. The plasma is assumed to have weak stationary 
density inhomogeneity, and thermal motions are not considered. Different polarizations, 
the left and right circular polarizations, are degenerate in a homogeneous isotropic plasma 
medium Q. In the presence of inhomogeneity, we show that this double degeneracy is lifted 
by a gauge connection term in the wave equation, which gives rise to Berry effects in the 
dynamics of polarized rays in the post geometric optics approximation. The polarization 
plane of a plane polarized wave rotates as a result of the geometric Berry phase, which is the 
Rytov-Vladimirskii rotation 

as 

. This is, of course, in contrast to the well known Faraday 
rotation, which is a dynamical effect due to the interaction of light and the magnetic field in 
a medium. Berry phase is a non-integrable phase factor ari sing from the adiabatic transport 
of a system around a closed path in its parameter space [lOj. Geometrically, it originates 
from parallel transport in the presence of a gauge connection in the parameter space 11]. 
When an electromagnetic wave travels in a weakly inhomogeneous medium, the direction 
of the wave vector varies slowly so that the parameter space in this case corresponds to the 
momentum space. The Berry curvature associated with the gauge connection in momentum 

n n 

space causes the optical Hall (or Magnus) effect [12], [13[], according to which, waves of 
left/right circular polarization deflect oppositely to produce a spin (polarization) current 
directed across the direction of propagation. Such Berry effects are very typical of spin 
transport and have been derived repeated 



in different inhomogeneous media [5 



'. y for various particles, in particular, for photons 
14KL2|. 

The paper is organized as follows. In section II, we write the plasma wave equation in 
an operator form by introducing a 'Hamiltonian' operator and show that, in the process 
of its diagonalization, a gauge connection emerges in the momentum space because of the 
inhomogeneity. In section III, we present the post-geometric approximation which is suit- 
able for weak inhomogeneity, and focus on the circularly polarized states by projecting the 
Hamiltonian on the polarization subspace. Finally, in section IV, we derive the Berry effects 
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in the dynamics of polarized rays by considering the post-geometric optics Hamiltonian, and 
establish the Rytov and the optical Hall effects in the plasma medium. 



II. PLASMA WAVE EQUATION AND THE GAUGE CONNECTION 

Consider an inhomogeneous plasma which can be generated by gravity or a position 
dependent electric field. Generally, both of them can exist and the balance equation of the 
forces in steady state is 

Pao{x)q a E (x) + p a0 (x)F ga - m a Vp a0 = 0, 

where p a0 is unperturbed mass density, F g is the gravitational force and a — e, i represent 
electrons and ions. Now, we want to derive the equation which determines the wave proper- 
ties of an inhomogeneous unmagnetized plasma. The procedure is similar to the derivation 
of the wave equation in homogeneous plasma [3] but we should notice that the density p a o( x ) 
is now a function of x. Considering small harmonic perturbation about steady state 

p a (x) = p a o(x) + p al (x)e~ luJt , 

V a = V al (x)e-^\ 
E = E (x) + E 1 (x)e- iuit , 
B = Bi(a;)e~ iw *, 

we linearize the two fluid and Maxwell's equations. The linearized momentum equation is 

-lUp a o{x)V al {x) = E X {X) +Pal{X) — r- — 7 V 



m a n a0 {x) \ p a0 {x) 

In deriving the above, the balance equation and p a \ = 1PaoPai(x)e~ iUlt / p a o(x) have been 

used. After combining it with linearized Maxwell's equations, the wave equation will be 

_ _ „ ,2 A u pe( x ) W pi( a3 )\ r 4tt_ g Q [ ( p al {x)p a0 \ p al (x 
V x V X E x - k [ 1 — \E 1 = — — 



Pa0{x) J Pa0{X) 



where k = uj/c and Up a (x) = 4Trp aQ (x)q^. The cold plasma approximation and assuming 
stationary and inhomo gene ous density is usually used to describe wave propagation in in- 
homogeneous media j^.l2C)l]. Under the cold plasma approximation, we can drop the right 
hand side of the equation and it takes the form 



V xVxJl = IJ|l-^-^|4. (1) 



Introducing 



e[x,u) 



CO 2 



and the dimensionless 'momentum' operator 



2 =7i (aj,w), 

-i^o V, (1) reads 



HEi = 0, 



(2) 



where the 'Hamiltonian' H is a matrix-valued differential operator with elements 



Hij(x,p,uj) = [p -n (x,u)]5ij - PiPj. 

In inhomogeneous plasma media the refractive index is also a function of frequency as well 
as coordinate. The refractive index is real for frequencies approximately above the plasma 
frequency and is imaginary for frequencies below it which causes the first to remain and the 
later to be absorbed in the media. 

The momentum operator obeys the standard commutation relations 

which has an important consequence with regard to diagonalization of H. To this end, we 
build a unitary matrix R(p) from the eigenvectors of the non-diagonal part PiPj, according 
to 



R(p) 



Py 



PxPz 



Px \ 



_ Px 



\ 



y/Px+Py 





p\JpI+pI p 

PvPz Pv_ 

P y/pj+P l p 

V'Pl+V 2 y 2* 

P P 



(3) 



where p = \p\, of course. The unitary transformation H — > R^HR, E\ — > RjE\ applied to 
equation (2), then yields 

H(x,p,u) =p 2 l - A-R^n 2 R, (4) 

where A = diag(0, 0, p 2 ) and 1 stands for the unit matrix. Of course, if x and p were classical 
(commuting) variables, the 'potential' part of the Hamiltonian (5) would simply reduce to 
n 2 l, and H would be completely diagonalized. However, in view of their anti-commutation, 
we have for the potential term in the momentum representation, 



R- 1 (p)n 2 (ik 1 V p )R(p) = n 2 {lik^ l V p - k G l A) = n 2 (xl - k^A) 
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where 

iA(p) = RT l VpR. 
In deriving the above, we have made use of the identity [ij], 

G-\x)f(d x )G(x) = f(d x + G-'dxG), 
since R~ X R=1, the vector matrix A is Hermitain. Thus 

H(x,p,u) = p 2 l - A - n 2 (xl - k^A). 

The canonical (generalized) coordinate conjugate to p is x which corresponds to the usual 
derivative ik^ 1 V p . In the absence of inhomogeneity, this canonical coordinate is the physical 
(observable) coordinate. With inhomogeneity present, however, the new coordinate, xl — 
k$ l A, corresponds to ikQ 1 T> p , where T> p is the covariant derivative defined by 

V p = lV p + iA(p). 

The inhomogeneity can be, thus, viewed as inducing the non-Abelian gauge connection 
(potential) A(p) in the momentum space. Such a gauge potential is a well known feature 



of spin transport ay , Il6l . l2ll |22J. A is a pure gauge potential, i.e., the corresponding 
field strength (curvature), V p x A + iA x A, is identically zero. The covariant derivatives, 
thus, commute so that the new coordinates also commute. Furthermore, the latter satisfy 
the same commutation relations with p as the canonical coordinates x. 

We consider that the inhomogeneity is weak, thus the direction of the wave propagation 
varies slowly (infinitely slowly for adiabatic density variations). Let us take this direction to 
be along the z-axis so that E iz is small (negligible in the adiabatic limit) compared to the 
other two components. Since we are considering polarization transport, we introduce the 
unitary matrix 

(l 1 \ 

% -i 
v V2j 

and write the wave equation in the helicity basis via the unitary transformation 

H^V ] HV = p 2 l- A-n^xl-k^A), E 1 -»■ E 1 = (E 1+ , E x _, E lz ) (5) 

where A = V^AV, and E\ a = (Ei x —iaEi y ) / a/2 represents the two circularly polarized states 
with helicity a = ±1. The elements i,j = 1,2 of H, thus, correspond to the polarization 
subspace, on which we shall focus. 



III. THE POST-GEOMETRIC APPROXIMATION 



We now proceed to study the dynamics imposed by the Hamiltonian <^ based on the 
Taylor expansion of n 2 (xl — A; ' 1 A). serves as a small parameter for the expansion, 
provided it is smaller than the length scale, L, of the density variations in the plasma. The 
zeroth order approximation (the geometric approximation), k^ 1 <C L, thus holds appropriate 
for large L, i.e., for a practically homogeneous plasma. In the zeroth approximation the 
Hamiltonian ([5]) becomes diagonal, 

i? (0) (cc,p, to) = diag(p 2 - n 2 ,p 2 - n 2 , -n 2 ), 

which implies a double degeneracy in the polarization subspace: electromagnetic waves 
with left/right circular polarizations have the same dispersion in a homogeneous isotropic 
plasma [7|. Furthermore, E\ z = 0, and the wave travels unrefracted along the z-direction, 
as expected. 

In inhomogeneous plasmas, the refractive index gradient or equivalently the plasma fre- 
quency gradient removes the polarization degeneracy. In the first approximation (the post- 
geometric approximation), which is suitable for weak inhomogeneity, the Hamiltonian (JHJ) 
takes the form 

H {1 \x,p,u) = H {0) (x,p) + k^Vn 2 ■ A. 

Calculating A via Appendix A, the non-diagonal correction term is seen to couple E a to 
E\ z , too. Since E\ z is small (negligible in the adiabatic limit) this coupling can be ignored. 
Consequently, we can project on the polarization subspace with the result 

H {l \x,p,u) = (p 2 -n 2 )l + kQ l Vn 2 .A ± a 3 , 

where A±(p) = p ^ p i^_ p ^ (~P2,Pi, 0) and <r 3 is the Pauli matrix. The wave equation, thus, 
breaks down into two independent equations for the left/right circularly polarized waves 
according to %„E a = 0, where 

H„ = l -[{p 2 - n 2 (r, co)} = \{p 2 - n 2 + ak^Vn 2 .A ± ). (6) 

(The factor | has been introduced for later convenience.) As remarked, the refractive- index 
gradient is responsible for the removal of polarization degeneracy. A±a 3 , or equivalently 
its eigenvalue Aj_a, is the (Abelian) gauge potential that emerges in the momentum space 



in the adiabatic approximation. The corresponding field strength (Berry curvature) and 
the physical coordinates are, thus, given by V p x A±a = —p~ 3 pa and r = x — k^ 1 A±a 
respectively. Note that because of the non-vanishing gauge field strength, which is the field 
of a magnetic monopole of charge —a situated at the origin of momentum space, the physical 
coordinates, now, do not commute: 

['"tj fj] = i&ko ^ijk^- 

IV. BERRY EFFECTS IN THE DYNAMICS OF POLARIZED RAYS 

Having constructed the Hamiltonian, the equations of motion of an electromagnetic ray 



can be obtained via the Hamilton's equations 23[. In the post geometric approximation, 



the semi-classical equations of motion of a circularly polarized ray are, therefore, 

where dot denotes derivative with respect to the ray parameter s, defined in terms of the 
ray length, I, by dl = nds. The physical coordinate r = x — k^ 1 A±a and the momentum 
p = k^ 1 k are now considered classical, of course. Along the trajectory, they represent the 
ray's position and (dimensionless) wave vector, respectively. Thus, using (jSJ), 

p=\v r n\ r=p + ak l1 ^-. (7) 
z p 

These, of course, reduce to the standard ray equations of geometric optics in the zeroth 
('classical') approximation, k^ 1 — > 0, where the left/right circularly polarized rays follow 
the same trajectory. However, in the post-geometric (semi-classical) approximation, as seen 
from (171) , the rays split due to the effect of Berry curvature of the momentum space (the 
magnetic monopole-like gauge field strength). The deflections from their classical (geometric 
optic) trajectories are given by 

p x dp 



5r = ak 1 / 



c p 3 



where C is the ray trajectory in momentum space. The resulting displacements are, there- 
fore, opposite and locally orthogonal to the direction of propagation. This, which is a general 
feature of spin transport, constitutes the optical Hall effect in inhomogeneous plasmas. 
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The phase change suffered by the ray in the course of its propagation is given by 



ut 



J k ■ dx = cut — k J p ■ dr + a J A ± ■ dp, 

where in the last integral, we have used the fact that Aj_ ■ p = 0. This integral represents 
the geometric Berry phase, which is of opposite signs for the two polarizations. Therefore, 
the polarization plane of a plane polarized ray rotates through the angle 

J A ± - dp = J cos 9 dcp, (8) 

where 9(ip) is the zenith (azimuth) angle in the spherical polar coordinates of the momentum 
space. This is the Rytov-Vladimirskii rotation. 



V. SUMMERY AND CONCLUSION 



It was shown that if we consider the post geometric approximation in unmagnetized 
plasma media with stationary slowly varying density, an Abelian gauge field (Berry con- 
nection) will result in the ray Hamiltonian. Appearance of this gauge field leads to an 
additional displacement of the photon of distinct helicity in opposite directions normal to 
the ray. Also, for a linear polarization, the rotation of the polarization plane occurs in inho- 
mogeneous unmagnetized plasmas. This rotation differs from the faraday rotation which is 
due to the presence of magnetic field in plasma media. In fact, it is shown theoretically that 
the rotation of polarization plane occurs even in unmagnetized plasma media as a result of 
Berry topological phase. 



VI. ACKNOWLEDGEMENT 



One of the authors (R.T.) wishes to thank to Dr. A. Chakhmachi for helpful discussions. 

Appendix A 

Components of gauge connection 

Direct calculation via (3) yields the following expression for gauge connection components: 
/ o p y Pz - Pv \ ( fi PxPz Px \ 



%A-y = 



PyPz 
Py 



p{pI+pI) 




Pyjpl+P 

PxPz 



PxPz 



\ Py/pl+Py P 2 ^Pl+Pl 



P 2 ^Px+Pl 








PxPz 

p(pI+pI) 


Px 

p^/pI+pI 


PxPz 





PyPz 


P(Px+Py) 


p 2 ^/pI+pI 


Px 


PyPz 





Pyfpl+Pl 


p 2 ^pI+pI 
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